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We investigate transport through a single-level quantum dot coupled to noncollinearly magnetized 
ferromagnets in the presence of localized spins in either the tunnel barrier or on the quantum dot. 
For a large, anisotropic spin embedded in the tunnel barrier, our main focus is on the impurity 
excitations and the current-induced switching of the impurity that lead to characteristic features in 
the current. In particular, we show how the Coulomb interaction on the quantum dot can provide 
more information from tunnel spectroscopy of the impurity spin. In the case of a small spin on the 
quantum dot, we find that the frequency-dependent Fano factor can be used to study the nontrivial, 
coherent dynamics of the spins on the dot due to the interplay between exchange interaction and 
coupling to external and exchange magnetic fields. 
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I. INTRODUCTION 

Recently, there has been growing interest in spin- 
dependent transport through nanostructures due to pos- 
sible applications for spintronics devices. Transport 
through quantum dots coupled to ferromagnetic elec- 
trodes is particularly interesting due to the interplay 
of strong Coulomb interaction on the quantum dot and 
nonequilibrium physics. Such systems have been realized 
experimentally, e.g., by coupling self-assembled semi- 
conducting quantum dots, 1-5 small metallic grains, 6-10 
quantum dots denned in InAs nanowires, 11 carbon nan- 
otubes, 12-16 or even single Ceo molecules 17 to ferromag- 
netic electrodes. Quantum dots coupled to ferromagnets 
have been also studied extensively from a theoretical 
point of view. 18-31 Especially interesting are quantum- 
dot spin valves, i.e., a single-level quantum dot coupled 
to noncollinearly magnetized electrodes. 32-39 Here, spin- 
dependent tunneling leads to a spin accumulation on the 
dot that in turn influences the transport through the sys- 
tem by blocking the current. In addition, there is an ex- 
change field acting on the dot spin. It is due to quantum 
charge fluctuations on the quantum dot and relies on the 
strong Coulomb interaction on the dot. It gives rise to 
a precession of the accumulated spin, thereby lifting the 
spin blockade of the dot. The resulting interplay of spin 
accumulation and spin precession leads to a number of in- 
teresting features in the transport characteristics, e.g., a 
shift of the peaks in the linear conductance with the an- 
gle enclosed by the magnetizations, 32,33 a broad region 
of negative differential conductance in nonlinear trans- 
port 33 as well as to characteristic features in the finite- 
frequency Fano factor 35 and a splitting of the Kondo res- 
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onance. 

While in the above studies all the spin dynamics takes 
place in the singly occupied orbital of the quantum dot, 
more complex spin dynamics appears when additional 
spin excitations are possible. These may be spin waves 
in the ferromagnetic leads, as discussed in Ref. 40. In 
the present paper, we investigate a different situation, 
namely, the coupling to a magnetic impurity either with 



a large, anisotropic spin or with a spin 1/2. 

We consider two different scenarios. In the first sce- 
nario, a magnetic impurity with a large, anisotropic spin 
is embedded in one of the tunnel barriers of a quantum- 
dot spin valve. Here, our main focus is on the spec- 
troscopy of the impurity spin as well as on its switching 
by the spin-polarized current. We consider the impu- 
rity spin in one of the barriers as this leads to a sim- 
pler conductance pattern (only processes involving this 
particular barrier can excite the spin), thereby simplify- 
ing the analysis of the impurity spin behavior. Our sys- 
tem is somewhat related to the case of transport through 
single tunnel barriers containing a magnetic atom or a 
single molecular magnet, that has been investigated ex- 
tensively in the recent past, both from a theoretical 41-49 
as well as from an experimental 50-56 point of view. It 
was shown that the steps observed in the differential 
conductance can be used to extract magnetic properties 
such as anisotropies of the atomic spin. 44-47,53 Further- 
more, the influence of nonequilibrium spin occupations 
was discussed, 49 explaining the overshooting observed at 
the conductance steps in experiment and predicting a 
super-Poissonian current noise. The absence of certain 
nonequilibrium features in turn was interpreted in terms 
of an anisotropic relaxation channel. For systems with 
magnetic electrodes, the possibility to switch the embed- 
ded spin by the spin-polarized current through the barrier 
was predicted theoretically 42,43 ' 48 and observed in exper- 
iment. 56 

In the model studied in this paper, the tunnel barrier 
containing the magnetic impurity connects a lead with 
a quantum dot. The transport behavior is, in this case, 
more complex since the spin dynamics of the embedded 
impurity is coupled to the charge and spin degrees of free- 
dom of the quantum dot. This has several consequences. 
First, interference between direct and exchange tunnel- 
ing through the barrier plays already a role for nonmag- 
netic electrodes in transport to lowest order in the tunnel 
coupling. This contrasts with the simpler case of a sin- 
gle tunnel barrier with a magnetic impurity, for which 
this interference only contributes for ferromagnetic elec- 
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trodes or in higher order transport. ' 48,56 ' Second, we 
demonstrate that the Coulomb charging energy of the 
quantum dot can help to perform tunnel spectroscopy 
on the embedded spin. Even if the excitation energy 
between the ground state and the first excited state of 
the impurity spin is larger than any other spin excitation 
energy, all spin excitation energies are accessible when 
additional charge states of the dot contribute to trans- 
port through the system, which is not possible for the 
single-barrier case. Third, we discuss a current-induced 
switching of the impurity spin. The impurity state influ- 
ence the spin accumulation on the quantum dot, which 
in turn acts back on the current through the system, 
leading to current oscillations as a function of the ap- 
plied bias. Interestingly, these phenomena occur even for 
small polarizations of the leads. While the current is only 
sensitive to the average value of the spin, we find that 
the zero-frequency Fano factor also contains information 
about the spin dynamics. Finally, fourth, we point out 
how monitoring the exchange field peak in the frequency- 
dependent Fano factor can detect the switching of the 
impurity spin for noncollinear magnetizations. 

In the second scenario, we consider a S = 1/2 im- 
purity side coupled to the spin of the electron on the 
quantum dot. Here, our main aim is to describe the co- 
herent dynamics of the two spins on the quantum dot. 
We consider the case that the impurity spin is located 
on the quantum dot, as this allows us to take the ex- 
change coupling between electron and impurity spin into 
account exactly. This model can serve to describe differ- 
ent situations. First, it can describe transport through a 
quantum dot that is doped with a magnetic atom. Trans- 
port through a quantum dot doped with a single man- 
ganese atom has already been studied theoretically. It 
was shown how the frequency-dependent shot noise can 
reveal the spin relaxation times. 58 Furthermore, the elec- 
trical control of the manganese spin state as well as the 
back action of the spin state on transport have been in- 
vestigated in the absence of Coulomb interaction in the 
quantum dot. 59 

Second, our model can be used to describe the cou- 
pling of the electron spin on the dot to a nuclear spin via 
the hyperfine interaction. In general, such a coupling is 
disadvantageous as it leads to decoherence of the electron 
spin and therefore can lift, e.g., the Pauli spin blockade in 
a double quantum dot. 60-63 However, it can also be used 
to dynamically polarize the nuclear spins in the quantum 
dot which in turn may be used to control and manipulate 
the electron spin. 64-68 

Transport through a quantum dot with a side-coupled 
spin 1/2 was discussed in Ref. 69 for the case of non- 
magnetic electrodes. It was shown how to extract the 
system parameters such as the exchange couplings, the g 
factors and spin relaxation times from measurements of 
the current and Fano factor. 

The case of noncollinearly magnetized ferromagnetic 
electrodes was recently investigated by Baumgartel et 
al. 70 for a large ferromagnetic exchange interaction. It 



was shown that in addition to a spin dipole, a spin 
quadrupole moment accumulates on the quantum dot, 
driven by a quadrupole current. 

In this work, we focus on the opposite regime of small 
exchange interaction between the side-coupled impurity 
and electron spin. This situation is particularly suited 
for the description of the weak hyperfine interaction. We 
discuss how the frequency-dependent Fano factor can be 
used to experimentally access the strength of the ex- 
change coupling for large and small external magnetic 
fields. Furthermore, for the case of a weak external mag- 
netic field, we show how the exchange field acting on the 
electron spin (but not on the impurity spin) gives rise to 
a highly nontrivial spin dynamics that manifests itself in 
the frequency-dependent Fano factor. 

Our paper is organized as follows. In Sec. II, we present 
the models that describe a magnetic impurity with a 
large, anisotropic spin localized in one tunnel barrier and 
a small spin localized on the dot, respectively. We intro- 
duce the real-time diagrammatic technique 71-74 that we 
use to calculate the transport properties in Sec. III. We 
discuss the form of the reduced density matrix of the 
quantum dot system and the generalized master equa- 
tion it obeys for the two systems under investigation in 
Sec. IV. Our results for the transport properties are pre- 
sented in Sec. V for an impurity in the barrier and in 
Sec. VI for an impurity on the dot. Finally, we conclude 
by giving a summary and comparing our results for the 
two models in Sec. VII. 



II. MODEL 

In this paper, we consider transport through a 
quantum-dot spin valve, i.e., a single- level quantum dot 
tunnel coupled to noncollinearly magnetized ferromag- 
netic electrodes. We consider additional magnetic impu- 
rities either with a large, anisotropic spin localized in the 
tunnel barrier or with a spin S = 1/2 on the quantum 
dot itself. In the following, we define the Hamiltonians 
for these two cases, respectively. 

A. Model A: Large spin in the barrier 

Our model A, which is schematically shown in Fig. 1, 
consists of a quantum-dot spin valve with an impurity 
embedded in the right tunnel barrier. As already men- 
tioned in Sec. I, we choose the impurity to be in the tun- 
nel barrier, as this will lead to a simpler conductance 
pattern, see Sec. V below, which allows to study the spin 
excitations more easily. In this case, the Hamiltonian can 
be written as the sum of four terms describing the two 
electrodes, the quantum dot, the spin and the tunneling 
between dot and lead, 

H = 2^ ^ r ^ dot -^spin + H tvi n,r- (2-1) 

r=L,R r=L,R 
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FIG. 1. (Color online) Schematic of model A, a quantum- 
dot spin valve with a magnetic impurity localized in the right 
tunnel barrier. 




FIG. 2. (Color online) Schematic of model B which consists 
of a quantum-dot spin valve with a magnetic impurity on the 
quantum dot. 



The first term, 



coupling to the left lead, it takes the form 



//,. 



k 



(2.2) 
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describes the ferromagnetic electrodes in terms of nonin- 
teracting electrons at chemical potential fi r . We quantize 
the electron spin in the direction of the magnetization 
of the respective lead. The spin polarization is defined 
as p r = (p r + — Pr-)/{pr+ + Pr~), where p rr7 is the con- 
stant density of states for majority (a = +) and minority 
(a = — ) spin electrons. 

The quantum dot is described by 



H, 



dot 



eel c a + Uc^c^ 



(2.3) 



Here, the first term characterizes the single, spin- 
degenerate quantum dot level with energy e measured 
with respect to the Fermi energy of the leads in equilib- 
rium. The charging energy U is needed to occupy the 
quantum dot with two electrons. As indicated in Fig. 1, 
we choose the quantization axis of the dot parallel to the 
magnetization of the right electrode as this simplifies the 
expressions for the tunnel Hamiltonian, see below. 

The third term in Eq. (2.1) describes the magnetic im- 
purity embedded in the right tunnel barrier as a localized 
spin with Hamiltonian 



-DSt~BS z . 



(2.4) 



We model the spin of magnitude S as having a uniaxial 
anisotropy D which favors the spin to be in the eigen- 
states pointing along the z axis. We, furthermore, as- 
sume the presence of a magnetic field B acting on the 
impurity spin. For concreteness, we assume this field to 
be pointing along the z direction. This choice is moti- 
vated by the presence of stray fields from the ferromag- 
netic electrode which have the tendency to align the im- 
purity along the magnetization of the electrode. As for 
the quantum dot, we quantize the impurity spin along 
the direction of the magnetization of the right electrode. 

The last part of the Hamiltonian (2.1) describes the 
tunneling between the dot and the electrodes. For the 
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LkH 



cos-c T -sm-c ; 
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sin -c T + cos -ci 



h.c, (2.5) 



i.e., majority and minority spin electrons of the leads 
couple to spin up and spin down states of the quantum 
dot due to our choice of quantization axes. The coupling 
to the right lead consists of two terms, 

-fftun.R = 2J i R 8 Rka C " + -?R a RkQ S " CafiCp + h.C. 
ktr 

(2.6) 

Here, cr a p denotes the vector of Pauli matrices. The first 
part describes direct tunneling between the dot and the 
leads. The second term describes exchange scattering 
from the impurity spin. 

The tunnel matrix elements and % (which can 
be chosen real) are related to the spin-dependent tun- 
neling rates via T rtT = 2ir\t r \ 2 p a . The total tunnel 
coupling is then given by T r = '^2 t7 F ra /2. Similarly, 
for the exchange tunneling, we relate the correspond- 
ing tunneling rate to its (real) matrix element by Jr = 
^Ij'rJ 2 P+ % P ~ ■ Furthermore, there will be a contribution 
due to the interference between direct and exchange tun- 
neling through the right barrier. It is characterized by 
i]\/TiiJft. Here, rj = ±1 determines the sign of the in- 
terference contribution which is governed by the relative 
sign of £r and jr. The upper sign, rj = +1 applies for 
equal signs of £r and j'r while the lower sign, r\ = — 1 
applies for different signs of % and jr. 



B. Model B: Small spin on the dot 

Model B consists of a quantum-dot spin valve with an 
additional spin localized on the quantum dot as shown 
schematically in Fig. 2. Here, we restrict ourselves to the 
case of a S = 1/2 impurity spin for two reasons. First 
of all, this keeps the size of the Hilbert space small while 
still giving rise to the nontrivial spin dynamics we are 
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interested in. Second, for a larger spin that additionally 
has some anisotropy, spin states would not be degenerate 
any longer, thereby destroying the possibility to observe 
the coherent spin dynamics, see also the discussion in 
Sec. III. The total Hamiltonian now takes the form 



H= H r + H, 



dot 



(2.7) 



=L,R 



describing the electrodes, the dot containing the impurity 
spin and the tunneling between dot and leads. The first 
part, H r , is identical to the one given in Eq. (2.2). For 
the second term, we have 



H dot = Y ec ^ c ° + ^ ( n t ~ n l) + Un t n i 



+ BS z + jJ2ci S ^c a ,, (2. 



The first two terms describe the single dot level with en- 
ergy e and Zceman splitting B due to an external mag- 
netic field. For simplicity, we assume the magnetic field 
to point along the quantization axis of the dot which 
we choose perpendicular to the magnetizations of the 
leads. The third term describes the Coulomb interaction 
U which is needed to doubly occupy the quantum dot. 
The fourth term describes the Zeeman energy of the ad- 
ditional spin on the quantum dot. Here, we assume the 
same g factor for the electrons on the dot and the im- 
purity spin, for a discussion of a system where electron 
and impurity spin have different g factors, cf. Ref. 69. 
Finally, the last term describes an exchange interaction 
between the spin of the electron on the dot and the im- 
purity spin. 

The eigenstates of the dot Hamiltonian (2.8) and their 
corresponding energies are summarized in Table I. The 
eight states can be classified according to the number of 
electrons on the dot. For the empty and doubly occupied 
dot, we have two states each that differ in energy by the 
Zeeman energy B with the impurity spin pointing up 
or down. For the singly occupied dot, we have in total 
four states, three triplet (S — 1) and one singlet state 
(S = 0). While the triplet states are energetically split 
by the Zeeman energy, singlet and triplet are split by the 
exchange coupling J. 

The coupling between dot and leads is described by 
the tunneling Hamiltonian 



rk 



*rk- 



cA +h.c, (2.9) 



where 0l = — 4>r = <j>/2 denotes half the angle enclosed 
by the magnetizations. 33 We relate the tunnel matrix ele- 
ments t r to the tunnel couplings T riJ as for the model dis- 
cussed in Sec. II A. Instead of using the tunnel coupling 
strength to the left and right lead, we can characterize the 
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TABLE I. Eigenstates of the dot Hamiltonian (2.8) and cor- 
responding eigenenergies. The first entry in each ket denotes 
the state of the quantum dot while the second entry charac- 
terizes the impurity spin state. 



dot-lead coupling alternatively by the total tunnel cou- 
pling r = Tl + Tr, and the asymmetry a = (IY — Tr) /T 
with — 1 < a < 1. 



III. TECHNIQUE 

In order to evaluate the transport properties of the two 
systems under investigation, we make use of a real-time 
diagrammatic technique 71-74 and its extension to sys- 
tems with noncollinearly magnetized ferromagnetic elec- 
trodes. 32,33 The basic idea of this approach is to inte- 
grate out the noninteracting degrees of freedom of the 
electrodes. We, then, arrive at an effective description of 
the remaining, strongly interacting subsystem in terms 
of its reduced density matrix p rcd . 

We denote by P^ 1 = (xi|p red |X2) the elements of the 
reduced density matrix, where |xi) and |xa) are eigen- 
states of the reduced system. The time evolution of the 
reduced density matrix elements is governed by a set of 
generalized master equations, 



e )P Xl 



X'xX'2 



yrrXiXipXi 
X2X2 X2 ' 



(3.1) 



The first term on the right-hand side describes the co- 
herent evolution of the reduced system. The second term 
describes the dissipative coupling to the electrodes. The 

generalized transition rates W XlX f are defined as irre- 

to X2X2 

ducible self-energies of the dot propagator on the Keldysh 
contour. They can be evaluated diagrammatically in a 
perturbation expansion in the tunnel coupling strength 
r. The corresponding diagrammatic rules are summa- 
rized in Appendix A. 

Expanding the density matrix elements as well as the 
generalized transition rates in a power series in the tunnel 
coupling r, we find that in the stationary limit the master 
equation for the off-diagonal matrix elements to leading 
order in the tunnel couplings takes the form 







< £ x 



)pXl 

il 1 X2 



(3.2) 
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if e Xl — e X2 >• T. As a consequence, coherent super- 
positions between states whose energy difference is large 
compared to the tunnel coupling have to be neglected 
in the sequential tunneling regime. On the other hand, 
for superpositions that satisfy e Xl 
equation in the stationary limit takes the form 



e X2 < r, the master 



= -i( E> 



£ )P Xl 



V w xix '} 

Z-^i X2X2 



pXi 

x'2 



(3.3) 



Xix 2 



to lowest order in the tunnel coupling. Here, the general- 



ized transition rates have to be evaluated at e Xl — s X2 







in order to consistently neglect all effects of order T 2 . 
Hence, in this case, the coherences will not vanish in 
general. 

We define the current through the system as the av- 
erage of the currents through the left and right tunnel 
barrier, / = (7 L — Ir)/2. It is given by 



I = ^e T W J P. 

2h 



(3.4) 



Here, we introduced the vector P which contains all den- 
sity matrix elements written as a vector to allow for a 
compact notation. The vector e T is a vector containing 
a 1 if the corresponding entry in P is a diagonal density 
matrix element and a otherwise. Finally, the quan- 
tity Wf contains the current rates W IxiXl which can be 
obtained diagrammatically similarly to the generalized 

transition rates W XlXl by replacing one tunneling vertex 

X2X2 

by a current vertex. The corresponding diagrammatic 
rules are given in Appendix A. 

The frequency-dependent current noise is defined as 
the Fourier transform of the symmetrized current-current 
correlation function S = (J(t)7(0) + J(0)7(t)) - 2(I) 2 . In 
the sequential tunneling regime and for low frequencies, 
^ S, r, we can write it as 



w 11 + w 1 (iioV) - w ) lw/ 

- 27r8(oj)(I} 2 + (uj^ -oj), 



P 

(3.5) 



IV. REDUCED DENSITY MATRIX AND 
MASTER EQUATION 

A. Model A: Large spin in the barrier 

For a quantum-dot spin valve with a large, anisotropic 
impurity spin embedded in the tunnel barrier, the eigen- 
states of the reduced system consisting of quantum 
dot and impurity spin are products of dot eigenstates 
|x) G {|0),| t)i| l)i \d}} and impurity spin eigenstates 
\S Z ) e{\ + S), I - S)}, \0 = \ X ) ® \S Z ). Assuming the 
energies of states with different impurity states to differ 
more than the tunnel coupling, Es z — 3> T, we have 
to neglect coherent superpositions between states with 
different impurity states. 

The reduced density matrix therefore takes a block di- 
agonal form given by 



/ Po, 











„rcd 
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m=-S 



\ 



P 



pi 

f,OT 





(4.1) 



P d , m ) 



In order to give a physically intuitive interpretation of the 
generalized master equations, we introduce the probabil- 
ities to find the dot empty, Po, m , singly occupied, P\. m — 
Pf,m + P\.,mi an d doubly occupied, Pd, m , with the impu- 
rity in state \m). We collect these quantities in the vector 
P,„ = (Po,m, Pi, mi Pd,m) T ■ We furthermore introduce 

the average spin on the quantum dot s x ,. 

1 and s, m = 



\<y.m ' : - '~~z.ui 2 ■ The set of 

master equations can then be split into one determining 
the occupation probabilities and one set governing the 
average dot spin. In the following, we will keep the time 
derivative on the left-hand side of the master equations 
explicitly to allow for a physically intuitive interpretation 
of the master equations. For the numerical discussion 
below, these derivatives are equal to zero, however. The 
master equations for the occupations are given by 



where W 7/ is obtained from W by replacing two tunnel 
vertices by current vertices. The frequency-dependent 
free dot propagator on the Keldysh contour is given by 



n (^ 



e X2 - £xi - w + i0+ 



(3.6) 



We stress that the frequency-dependent current noise for 
lj < r is only sensitive to coherences between states with 
e xi — £ X2 S r, i.e., we can savely neglect all other coher- 
ences in the calculation of S^w) as in the evaluation of 
the master equation and the stationary current. While 
our discussion of the real-time diagrammatic technique 
so far was rather general, in the next section, we turn 
to the explicit form of the density matrix as well as the 
generalized master equations for the two systems under 
investigation. 



W (0) P 

f w (0) p 



V {0) s 



n L 



W( +1 )p , , 4. w( _1 )p 
PI ^m+l 1 vv r r m— 1 



V (0) s 



Here, W, is a matrix which describes processes that 
transfer an electron between the dot and lead r and 
change the state of the impurity spin from |m) to \m+a). 
Since the impurity is located in the right tunnel barrier, 
tunneling through the left lead cannot change the impu- 
rity state and therefore w[ ±l ' = 0. Changing the im- 
purity state is possible, however, for tunneling through 
the right barrier which provides a coupling between P m 
and P m ±i- Similarly, is a vector which describes 

the coupling of the occupation probabilities to the spin 



n R + Ki +1) s m+ i • n R + V I 



(-1) 



Sra-l • nR. 

(4.2) 



G 



on the dot, a feature characteristic of a quantum-dot spin 

valve. Again, we have = 0. The precise form of 

W { r a) and Vr (a) is given in Appendix B 1. 

The time evolution of the dot spin obeys a Bloch-type 
equation, 



~df 
~dJ 



(0) 
acc,L 

(-i) 

acc.R 





(0) 




acc . R 




v (0) 


V dt , 


' rcl,L 




X (B 



~dt 
~dt 



(+i) 

acc,R 
(0) 

rel,R 



The hrst four terms on the right-hand side describe the 
non-equilibrium spin accumulation on the quantum dot 
due to the tunneling from and to the spin-polarized leads. 
Similarly to the master equation for the occupations, we 
get accumulation terms that change the state of the im- 
purity when tunneling takes place through the right bar- 
rier. 

The next two terms account for the relaxation of the 
spin on the dot due to the tunneling out of an electron or 



the tunneling in of an electron forming a spin singlet with 
the electron already present on the dot. As these terms 
arise from generalized transition rates which start and 
end in a singly occupied state, in the sequential tunneling 
approximation the state of the impurity spin cannot be 
changed in these processes. We give the explicit forms of 
the accumulation and relaxation terms in Appendix B 1. 

The last term describes the precession of the dot spin 
in an exchange field due to virtual tunneling to the leads. 
For the coupling to the left lead, we find the usual ex- 
change field 33 which is independent of the state of the 
impurity spin. It is given by 

B m . L = -n L ^($ L ( £ )-$ L ( £ + [/)), (4.4) 

7T 



where § r (x) = Rc^ f| 



+ i 



0{x- 



2tt 



and denotes the 



digamma function. While the first term arises from the 
spin-dependent level renormalization of an electron vir- 
tually tunneling to the lead and back, the second term 
stems from processes where an electron first tunnels onto 
the dot and then back into the lead. 

The exchange field due to the coupling to the right lead 
which points in the direction of n R is given by 



B 



m,R 



-n R ^ $ R (e) - $ R (e + U) 

7T 



+ n R 



A+(m-l)J R 



^*k( £ - a_) + ^$ R ( £ + U + a_) - In ^ 



n R - 



A_(m + 1)J R 



X — ^$ R (e + a+) + 4^* R (e + U - a+) - In 



~2V 



(4.5) 



where 7 = prT r + m 2 p R J R + 2mr}y/T R J R , a± = B + 
(2m ±1).D and A±(m) = S(S + 1) - m(m± 1). 

Here, the first term on the right-hand side describes 
exchange field contributions due to virtual tunneling be- 
tween dot and lead that does not change the state of 
the impurity spin. The other two terms are due to vir- 
tual tunneling where the intermediate state has an in- 
creased/decreased impurity spin state. These processes 
give rise to a logarithmic divergency of the exchange field 
cut off by the bandwith W of the lead electrons. To 
understand this, we consider the impurity in the state 
\S X = S). In this case, only the sequences | l,S) —> 
10, S - 1) -> I I, S) and I I, S) -> \d, S - 1) -> | |, S) of 
(virtual) transitions are possible while there are no such 
processes starting from | t;^)- Hence, these processes 
only renormalize the energy of the spin down state, giv- 
ing rise to the logarithmic divergency. Similarly, when 
the impurity is in any other state, the logarithmic con- 
tributions to the exchange field do not cancel between 
processes that increase and decrease the intermediate 



impurity spin state. As we only consider sequential 
tunneling, it is clear that our results are only valid if 
J R ln^ « feT and J R In ^ < r r , otherwise higher 
order logarithmic corrections become important. 



We, therefore, find that the presence of the impurity 
spin in the tunnel barrier has two basic effects on the 
exchange field. First of all, it alters its strength. Sec- 
ond, due to the presence of the spin-flip processes, it also 
alters its energy dependence (cf. Fig. 3), giving rise to 
additional peaks and dips whose separation is governed 
by the anisotropy D, the Zeeman energy B, and the size 
of the impurity spin S. Since the transition energies be- 
tween the various impurity states depend on the states 
itself, the position of the additional peaks and dips de- 
pends on the value of S z . 
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FIG. 3. Exchange field due to virtual tunneling through the 
right barrier as a function of bias voltage V for p = 0.3 (upper 
panel) and p = 0.9 (lower panel). Other parameters are e = 
-300fc B T, U = 250fc B T, D = 17.5fc B T, B = Wk B T, T R = 
IOJr, n = +l,W= 500fc B T, 5 = 3. 
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Parameters 


Superpositions of 


(i) 


b, j, \b - j\ > r 




(ii) 


s>r, j<r 


\T°), \S) 


(iii) 


s,j>r, |b- j\ <r 


IT'}, \S) 


(iv) 


j > r, b < r 


lot), |<H) 
|t+), |t°), |r-) 

|dt), |<U) 


(v) 


b, j < r 


lot), |0 4-> 

|T+), |T°), |T"), |5> 

Mt), |d;> 



TABLE II. Coherent superpositions that have to be taken into 
account in the sequential tunneling regime for different values 
of the external magnetic field B and the exchange coupling J 
between the two spins. 
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i.e., apart from the eight diagonal matrix elements that 
describe the probability to find the system in one of its 
eigenstates there can be up to 16 coherences. Coher- 
ences between states with different electron numbers van- 
ish due to the conservation of total particle number. As 
discussed above, depending on the energy differences of 
the states forming the coherent superposition, we either 
have to take them into account or neglect them in the 
sequential-tunneling regime. In Table II, we summarize 
the different transport regimes that arise consequently. 

In the following, we will only consider the cases (ii) 
and (v), i.e., we only consider the case of weak exchange 
couplings, J < T. When a large magnetic field is ap- 
plied, B>r, only S-T° coherences have to be taken into 
account. When the externally applied magnetic field is 
weak, B < T, we have to take into account all coherences. 
There are two reasons focusing on the two cases. On the 
one hand, they are particularly suited to demonstrate 
the information about the transport processes contained 
in the finite-frequency noise. On the other hand, the 
cases of small exchange couplings are suited to describe 
the influence of nuclear spins, that couple to the electron 
spin via hyperfine interaction, on transport through the 
quantum dot. 



1. Case (ii): Large magnetic field 

We first turn to the discussion of the master equation 
in case (ii) where B ^> T, J < T. In this case, there are 
only superpositions of S and T° present. It is therefore 
natural to introduce the isospin I via 



It = 



PrpQ H~ Pq 



P T o - P s 



pS _ pi 



Iy = f 



I Q = P T o+Ps, 



B. Model B: Small spin on the dot 



For the case of a S = 1/2 impurity localized on the 
quantum dot, the reduced density matrix in the most 



to bring the master equation into a physically in- 
tuitive form. Similar to the case of an ordinary 
quantum-dot spin valve, 33 we can now split the mas- 
ter equation into one set governing the occupation 
probabilities that we summarize in the vector P = 
{Pot,Poi,PT+,Io,PT-,Pdt,Pdi) and one set governing 
the time evolution of the isospin I = (J x , I y , I z ). How- 
ever, there is an important difference. While in the or- 
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dinary quantum-dot spin valve, there is a real spin accu- 
mulating on the quantum dot, here we have an isospin 
accumulation as a real spin accumulation is suppressed 
by the large external magnetic held. The master equation 
of the occupation probabilities reads 

P = W-P + V(I-0, (4.7) 

where W denotes a matrix that contains the golden rule 
transition rates between the various dot states and V is 
a vector that characterizes the influence of the isospin on 



the dot occupation whose precise form is given in Ap- 
pendix B 2. 

The master equation that governs the time evolution 
of the isospin is given by 

i =(£L+(§L+ i *? b - <"> 

Here, the first term, 



dt. 



= \ J2 r r l-f+(e - B/2)P^ + /+(e + B/2)P ai -f+(e + U + B/2)P dt + f-(e + U- B/2)P di 



f-(e - B/2) - f-(e + B/2) + f+(e + U + B/2) -f+( £ + U- B/2) 



Io}e x , (4.9) 



describes the accumulation of the isospin along the x axis 
due to electrons tunneling onto and off the dot. Similarly, 
the second term describes a relaxation of the isospin 



f~{e-B/2) + /-(e + B/2) 



relax 

+f+(e + U- B/2) + f+(e + U + B/2) 



I. (4.10) 



The exchange field describes the level splitting between 
|T°) and IS") which is due to the finite exchange coupling 
J as well as due to virtual tunneling processes that renor- 
malize the energies of the two states in a different way. 
As can be inferred from Eq. (4.8), it gives rise to a pre- 
cession of the accumulated isospin around the exchange 
field. 



Finally, the last term describes the precession of the 
isospin in an exchange field that is given by 



B rx = $ r (e - B/2) - $,,(e + B/2) 
Zir 



$ r (e + U + B/2) - $ r (e + U - B/2) 



B ry — 



(4.11) 



The current through the quantum dot is given by 



/ = 



/+( £ + B/2) + f+(e - B/2) )(Pot + Poi)-( + U + B/2) +f-( £ + U- B/2) )(P df + P di 
/ L " ( £ + B/2) -/+(£+ T B/2) ) /'•, -(./,>- B/2)- /-/(: + !' + B/2)) P r 
' ' / L -( £ + B/2) + f-(e - B/2) -f+(e + U + B/2) -f+(e + U- B/2) ) J 



(L-)R). (4.12) 



/ L - (e + B/2) - / L " (e - B/2) + f+( £ + U + B/2) -f+(e + U- B/2) I 



It depends on the occupation probabilities as well as on the x component of the accumulated isospin. This re- 
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sembles the normal quantum-dot spin valve where the 
current also depends on both, the dot occupations as 
well as on the dot spin. 33 

2. Case (v): Small magnetic field 

We now turn to the discussion of the master equation 
in the case B, J <T. In this case, we have to include all 
coherences of the reduced density matrix, Eq. (4.6). 

To allow for a physical interpretation of the different 
matrix elements, we introduce the probabilities to find 
the quantum dot empty, Pq, singly occupied, Pi, and 
doubly occupied, Pa- Furthermore, we introduce the ex- 
pectation values of the electron spin on the dot, Si, as 
well as the expectation values for the impurity spin when 
the dot is empty, So, singly occupied, S2, and doubly oc- 
cupied, Sd- While for a single spin 1/2 the description of 
its density matrix in terms of spin expectation values is 
sufficient, this is in general no longer true for a system 
of two spin 1/2 particles. 70 For the case of small mag- 
netic fields that we consider here, we therefore have to 



introduce in addition the expectation values of the scalar 
product between electron and impurity spin, Si • S2, and 
their vector product, Si x S 2 . Finally, we also need to 
introduce the quadrupole moment 70 ' 75 

Qij = 2 (SliS2j + SijS2i) — -Si • S2<%. (4.13) 

The quadrupole moment is a symmetric tensor, Qij = 
Qji. Its diagonal elements are not independent of each 
other as they satisfy the sum rule Y^iQa = 0' i- e -' tne 
trace of Q vanishes. 

In Appendix C we give the explicit expressions that 
relate the above quantities to the density matrix ele- 
ments in Eq. (4.G). We note that in the case B T, 
J < r where we only have taken into account the S-T° 
superpositions, we could have expressed the reduced den- 
sity matrix in terms of the quantities just introduced as 
well. However, by choosing a description in terms of the 
isospin, we obtain a much simpler master equation. 

Using the physical quantities we just discussed, we can 
split the master equation into several sets. The first set 



d_ 
db 



(Po 
Pi 
\P d 



-2f+(e) /-(e) 

2/+(e) -f-( e )-f+(e + U) 2f-(e + U) 
f r + (e + U) -2f~(e + U) 




2p r r r 



fr(e) 

-f-(e) + f+(e + U) I S r n r , 
-f+(s + U) 

(4.14) 



describes the evolution of the occupation probabilities. It 
takes a form identical to the case of a normal quantum- 
dot spin valve, i.e., it exhibits a coupling of the occu- 
pations to the spin accumulated on the quantum dot. 
Interestingly, the occupations do not couple neither to 



the impurity spin, the scalar or vector product of Si and 
S2 nor to the quadrupole moments directly. They are 
only influenced by these quantities due to their influence 
on the accumulated electron spin Si. 

The equation governing the time evolution of the elec- 
tron spin in the dot is given by 



rfSi 



E 



/v r, ( /+( £ )P - /r(£) fr + ( £ + u ) Pl _ f - {£ + U)P \ n? , _ ?! _ Sl x Br 



where l/r r = r r (/- (e) + /+ (e + U)) and B r , ox = 
— PEjl ($ r (g;) — $ r (e + U)) n r is the usual exchange field 
acting on the electron spin accumulated on the dot. 
Again, we find a strong similarity to the case of the nor- 
mal quantum-dot spin valve. While the first term in 
brackets describes the accumulation of spin on the dot 
along n r due to spin-dependent tunneling of electrons 
between dot and leads, the second term describes a re- 



-Si x B + J(Si x S 2 ), 

(4.15) 
I 

laxation of the dot spin due to tunneling. The third term 
in brackets describes the precession of the dot spin in the 
exchange field generated by virtual tunneling between 
dot and leads. The last two terms finally describe the 
influence of an external magnetic field and the exchange 
coupling to the impurity spin. 

The master equations for the impurity spin in the pres- 
ence of zero, one and two electrons on the dot can be 
written as 
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d_ 
dt 



( S 
S 2 



-2f+(e) /-(e) 

^Tr r 2/+(e) ~f-(e)-f+(e + U) 2f~(e + U) 

V fHe + U) -2f-(e + U) 

^2p r r r | -f-(s) + f+(s + U) 




Q+-Si-S 2 ) n 



r+^(SlX S 2 ) X 






/S \ 




) 


S 2 x B + J 


S 


ixS 2 


\s d ) 




/ 



(4.16) 



Here, the first term on the right-hand side describes tran- 
sitions between the three quantities by tunneling of elec- 
trons in analogy to the first term in the equation for the 
occupations, Eq. (4.14). The second term characterizes 
the coupling to the quadrupole moments as well as the 
scalar and vector product of Si and S 2 . This resembles 
the coupling of the dot occupations to the electron spin 



on the dot in Eq. (4.14). Finally, the terms in the third 
line describe the precession of the impurity spin in an ex- 
ternally applied magnetic field as well as the influence of 
the exchange interaction between electron and impurity 
spin. 

The master equations governing the time evolution of 
the scalar and vector product between the electron and 
impurity spin are given by 



dt 



(S lX S 2 )=]T 



Pr T r f+(e)8 a 



f-(e)-f+(e + U) 



Sa- f-(e + U)S d ) xn. 



Si x S 2 



Q--(Si-Sa) ) B r 



(SixS 2 )xB + -(Si-S 2 ), (4.17) 



!<s,.*)-x; 



Pr T r /+(e)S 



/-( £ ) + f+(e + U) 



S 2 - f-(s + U)S c 



n r — 



Si-S 2 



(Si x S 2 )-B r 



(4.18) 



Their structure closely resembles Eq. (4.15) in that there 
are terms which describe the accumulation, relaxation, 
and the influence of the spin precession due to the ex- 
change field. Furthermore, the vector product turns out 



to be sensitive to an external magnetic field as well as to 
the exchange coupling between the spins. 

Finally, the master equation for the quadrupole mo- 
ment takes the form 



Pr^rfr( £ ) Q i^O.i^j + ^0,]^-r,i) ~ ^(S • n r )Sij 

f-(e)-f+(e + U) (1 1 



Oil _ 1 n 
r r 2 I 2 



o ( x S 2 )i(B T . jex ) :) - + -(Si x S 2 ) J (B I , iCX ) i — -(Si x S 2 ) ■ B rjex <$y 



1 ,„ „ , , 1 



' ^iimQ/jBm £jlm. Qii-^m • (4. 19) 
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The first three terms on the right-hand side describe 
the accumulation of quadrupole moment on the quan- 
tum dot. Similarly, the fourth term is related to the 
relaxation of the quadrupole moment. Finally, the other 



terms describe the precesional motion of the quadrupole 
moment in the exchange field as well as due to an external 
magnetic field. 

The current through tunnel barrier r is given by 



I r = T r f+(e)P -T, 



f-(e)-f+(e + U) 



Pi-r r /-( £ + c/)p rf -pr r 



f-(e) + f+(e + U) 



Si ■ n r . 



(4.20) 



Although this is precisely the same form as for the normal 
quantum-dot spin valve, the current nevertheless con- 
tains information about the nontrivial spin dynamics on 
the dot, as the master equation for the dot spin couples 
to the other density matrix elements. 



V. RESULTS - LARGE SPIN IN THE BARRIER 

In this section we discuss the transport properties of a 
quantum-dot spin valve with a large, anisotropic impu- 
rity spin located in the right tunnel barrier. We will focus 
our attention on systems that are symmetric in the sense 
that pl = Pr = P and Tl = r R = T/2. Furthermore 
we assume the bias voltage to be applied symmetrically, 
V h = -V R = -V/2. 



A. Collinear magnetizations 

In the following, we are going to discuss transport for 
collinear magnetizations. Most of the time we will re- 
strict ourselves to the case of parallel magnetizations as 
we focus on moderate polarizations where the effects due 
to the impurity spin are much more important than the 
effects due to the relative orientation of the leads. 



1. Interference between direct and exchange tunneling 

From the form of the tunneling Hamiltonian (2.6) it is 
obvious that interference can take place between elec- 
trons tunneling directly through the barrier and elec- 
trons experiencing an exchange interaction with the im- 
purity spin. For transport through a single tunnel barrier 
containing a localized spin 41,43 ' 49 ' 76 ' 77 , the interference 
contributions cancel between the spin-up and spin-down 
channel. Only for ferromagnetic leads 48 ' 56,57 ' 78 or in the 
presence of spin-orbit interactions, one is sensitive to the 
interference terms. 

This is different for the system under investigation 
here. We find that the interference terms influence the 
current even for unpolarized leads. In contrast to the 
single-barrier case where we just have to sum up the con- 
tributions from spin-up and spin-down electrons to the 
current, in the quantum dot case we have to separately 



r 



compare the rates for tuneling in and out of the dot for 
each spin direction. While in the nonmagnetic case equal 
amounts of spin-up and spin-down electrons enter the dot 
from the left lead, the rates for leaving the dot are dif- 
ferent due to the interference terms. This in turn gives 
rise to a spin accumulation on the quantum dot which 
reduces the current through the quantum dot. 

Unfortunately, in our system there is no way to tune 
the phase of the interference terms experimentally as is 
possible, e.g., in an Aharonov-Bohm interferometer and 
thereby check the influence of the interference terms on 
the current. Nevertheless, it should be possible to detect 
the presence of the interference term experimentally and 
to detect its sign. Approximating the Fermi functions as 
step functions which is reasonable away from the thresh- 
old voltages, we can calculate the current through the 
system analytically in the various transport regions. 

We first consider the case of unpolarized leads, p r = 0. 
In this case, the current in region I where transport takes 
place through the singly and doubly occupied dot (cf. 
Fig. 4) is given by 



2r L (r R + ^ 2 j R ) 
r L + 2(r R + s 2 j R )' 



(5.1) 



i.e. it is sensitive to the couplings and the size of the 
barrier spin but not to the interference term. Similarly, 
in region II where spin excitations become possible, the 
current turns out to be insensitive to the interference 
term, 



2r L (r R + g(g + i)j R ) 
r L + 2(r R + s(s + i)j R )' 



(5.2) 



This is different in region III where transport takes place 
through the empty and singly-occupied dot but the spin 
cannot be excited. Here, the current is given by 



/in = 



2r L ((r R + s 2 j r ) 2 - 4S 2 VrW) 



(r R + s 2 j R )(2r L + r R + s 2 j r ) - 4S 2 yr R j R z ) 

(5.3) 

i.e., the current now also depends on the interference 
term. From Eq. (5.3), we infer that for T R = SJr the 
current vanishes exactly in region III. Equation (5.3) also 
shows that the current in region III is only sensitive to 
the absolute value of the interference term but not to its 
sign. 
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FIG. 4. (Color online) Schematic of the differential conduc- 
tance as a function of level position and applied bias voltage. 
Thick black lines mark the onset of transport through the dot. 
Blue dashed lines indicate the onset of impurity excitations. 
Dashed lines mark the gaining of energy from the impurity to 
allow all three charge states of the dot. 

This is different in the regime where the empty and 
singly occupied dot contribute to transport and spin ex- 
citations are possible. As the analytic result for the cur- 
rent in this regime is rather lengthy, we do not give it 
here. Instead, we now focus on transport for parallely 
magnetized leads. In region I, the current is now given 
by 

2T L (r R + S 2 J R + 2 P S V VT^h) 

1 r L + 2(r R + s 2 </ R ) ■ (& - > 

Here, the current is clearly sensitive to the sign of the 
interference term which provides a way to access it in ex- 
periments. Similar expressions for the current in regions 
II and III can be found for parallely magnetized leads. 
As these expressions are rather lengthy, we do not give 
them here. 



schematically in Fig. 4. For level positions e < — U , the 
dot is doubly occupied when no bias is applied. Upon in- 
creasing the bias, we enter region I where transport takes 
place through the singly and doubly occupied dot. When 
increasing the bias voltage above the blue dashed line, ex- 
citing the spin becomes possible, similarly as in a single 
barrier discussed above which only provides information 
about a linear combination of D and B. However, upon 
increasing the bias further, we reach the red dotted line. 
At this point, the electron on the dot with an energy be- 
low the right Fermi level can gain enough energy to leave 
the dot to the right lead by changing the impurity state 
from S z = S — 1 to S z = S. As this opens up a new 
transport channel, the onset of this process yields a sig- 
nal in the differential conductance. Similarly, at the next 
yellow line the process 1 .J.) ® | - S + 1) -> |0) ® | - S) be- 
comes possible, again giving rise to a conductance signal. 
This scheme continues for all transitions of the impu- 
rity spin which are characterized by transition energies 
(2S Z — 1)D + B for S z > and (-(2S 2 + l)D - B) for 
S z < 0. When we finally reach the thick black line, the 
empty dot state can also be reached by ordinary tun- 
neling events. For even larger bias voltages we find an- 
other series of conductance signals which are now associ- 
ated with electrons in the lower level which become able 
to excite the impurity when leaving to the right lead. 
The important difference to the small bias case discussed 
above is that now transport through the upper level can 
bring the impurity spin into all excited states such that 
transitions between these states also are all visible. 

For unpolarized leads, the conductance pattern dis- 
cussed above allows to determine D and the absolute 
value of B. As no spatial direction is distinguished, there 
is no possibility to determine also the sign of B. This is 
different for polarized leads where the spatial symmetry 
is broken by the magnetizations. In this case, one finds 
that the differential conductance shows an alternating 
pattern of positive and negative differential conductance 
which depends on the sign of B. We discuss the mecha- 
nism leading to this behavior in the following. 



2. Spin spectroscopy 

Inelastic spin tunneling spectroscopy 44-51 ' 53-55 ' 79 ' 80 al- 
lows to determine the spectrum of a spin embedded in a 
tunnel barrier by studying the steps in the differential 
conductance that occur whenever an inelastic transport 
channel opens up. However, for a simple spin Hamilto- 
nian of the form (2.4) this does not allow to determine the 
parameters D and B separately. In this case, the energy 
difference between the ground state and the first excited 
state is larger than all other excitation energies. 80 ' 81 In 
consequence, as soon as the system can be brought into 
the first excited state, all other excited states can also 
be reached energetically. Hence, there would be only a 
signal at A = E S -i - E s = (2S - l)D + B. 

This is different for the system additionally contain- 
ing a quantum dot between the electrodes, as is shown 



3. Current-induced switching and spin- dependent transport 

We now turn to our main result for the model contain- 
ing an impurity spin in the barrier, that is we discuss the 
interplay between the dot and impurity spin and its man- 
ifestation in the transport properties. In particular, we 
investigate how the the dynamics of the dot and impurity 
spin gives rise to the sequence of positive and negative 
differential conductance features at the red and yellow 
dotted lines in Fig. 4 for polarized electrodes and which 
is manifest in the current oscillations of the I — V char- 
acteristics shown in Fig. 5. (In our discussion, we focus 
on the case B > 0. For B < 0, the impurity ground state 
is | — S) and basically the role of the cases ?7 = +1 and 
rj = — 1 are interchanged.) 

As tunneling into a ferromagnet is spin dependent, 
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FIG. 5. (Color online) Current as a function of bias voltage for 
parallel magnetizations, p = 0.3 and Tl = Tr = Jr. Other 
parameters are the same as in Fig. 3. 
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FIG. 6. (Color online) Average z component of dot and im- 
purity spin as a function of bias voltage for rj = +1. Other 
parameters are the same as in Fig. 5. 



we find that the rates for the impurity spin transitions 
|m) — > \m + 1) and |m) — > \m — 1) are different, in gen- 
eral, because one transition involves tunneling of a mi- 
nority spin electron, while the other involves tunneling 
of a majority spin electron. We therefore have for the 
rates W m ^ m +i oc 1 — p and W m ^ m -i oc 1 +p. As a con- 
sequence, once exciting the impurity spin becomes ener- 
getically possible, the spin-polarized current through the 
right barrier has the tendency to flip the impurity spin 
into the state | — 5) which in Figs. 5 and 6 occurs at 
V = 2{e + U + A) = 295fc B T. 

At larger bias, transitions where energy is gained from 
the impurity come into play. Alternatingly, they cither 
lower the z component of the impurity spin, | — 5 + i) — > 
| - 5 + i - 1), or raise it, |5 - i) -> \S - i + 1), with 
i <E [0,5]. Hence, the expectation value of S z is found 
to oscillate as a function of bias voltage as is shown in 
Fig. 6. 

The other key ingredient for understanding the con- 
ductance oscillations is the fact that tunneling through 
the right barrier is spin dependent in two respects. On 



the one hand, there is the dependence on the spin of the 
tunneling electron due to the polarization of the lead, 
mentioned already above. However, the tunneling also 
depends on the state of the impurity spin as the tunnel- 
ing rate for spin-up (down) electrons is proportional to 
|£r + 5 2 jr,| 2 (|£r — 5 z jr,| 2 ). This kind of spin dependence 
is, then, responsible for relating the dot spin to the im- 
purity spin state. If, e.g., the impurity is in the ground 
state | + 5) , spin up electrons can leave more easily to 
the right lead than spin down electrons. This results in a 
spin down accumulating on the quantum dot, see Fig. 6, 
2(e + U) < V < 2(e + U + A). When the z component 
of the impurity spin has a negative expectation value, 
the situation is reversed. Now, spin down leaves the dot 
more easily such that spin up accumulates on the dot. 

The spin accumulation on the dot affects the current 
through the system. If the lower dot level is occupied by 
a spin up electron, the Pauli principle prevents a second 
spin up electron to enter the dot while transport of spin 
down electrons is suppressed by their smaller density of 
states. Hence, the accumulation of spin up on the dot 
suppresses the current. On the contrary, a spin down in 
the lower level does hardly affect the current, as it can 
proceed by spin up electrons tunneling through the upper 
level. 

Hence, in symmary we find that the interplay between 
the current-induced switching of the impurity spin and 
the spin-dependent tunneling through the right barrier 
results in an interesting spin dynamics in the quantum- 
dot spin valve which manifests itself in the transport 
properties of the system. We emphasize that the trans- 
port signatures discussed here are present for moderate 
polarizations. Hence, the polarizations of Fe, Co, or 
Ni (Refs. 17, 82) should be sufficient to experimentally 
detect them. For very large polarizations, the current 
oscillations discussed above are even absent, because in 
this case the current-induced switching mechanism is so 
strong that the impurity will be kept in state | — 5) once 
spin excitation becomes possible. 



4- Giant Fano factor 

In the above discussion of the current oscillations, we 
found the current to be sensitive to the average impurity 
spin only. The current noise however sheds more light 
on the dynamics of the impurity spin. As can be seen in 
Fig. 7, in the region where energy can be gained from the 
impurity to allow electrons to leave the dot to the right 
(corresponding to the region of dotted lines in Fig. 4), the 
Fano factor is strongly enhanced (note the logarithmic 
scale) . 

The mechanism which gives rise to these giant Fano 
factors is the following. As can be seen in the upper 
panel of Fig. 8, the probability to find the impurity in the 
states | ± 5) is finite for both states in the region where 
the Fano factor becomes large. The dot spin follows the 
behavior of the impurity spin which means that it will 
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FIG. 7. (Color online) Fano factor as function of bias voltage. 
Parameters are the same as in Fig. 6. The sharp peak at 
V = 295/cbT arises at the onset of spin excitations (dashed, 
blue line in Fig. 4). 
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FIG. 8. Probabilities to find the impurity in spin state \S Z ) as 
a function of bias voltage (upper panel) . The middle (lower) 
panel shows the probability to find the system in the state 
| t) ® (I i) ® l&z))- Parameters are the same as in Fig. 6. 



point either up or down, depending on the sign of S z as 
shown in the middle and lower panel of Fig. 8. The two 
dot spin configurations carry a different current, because 
only for the spin up accumulation transport through the 
dot becomes spin-blockaded. Hence, the system switches 
between two current states on a rather large time scale 
(as it takes several spin flip processes to reverse the im- 
purity spin) which gives rise to random telegraph noise. 
Similar behavior can be found, e.g., in arrays of move- 
able colloid particles, 83 transport through molecules with 
strong electron-phonon coupling and strong vibrational 



relaxation 84 or double dot Aharonov-Bohm interferome- 
ters. 85 

It is interesting to note that the giant Fano factor oc- 
curs for parallel as well as antiparallel magnetizations. 
This is due to the fact that the moderate polarizations 
chosen here the effects of bunching for parallel magneti- 
zations and spin blockade for antiparallel magnetizations 
are rather weak such that they are unimportant com- 
pared to the random telegraph switching. 

B. Noncollinear magnetizations 

So far, we only discussed the transport properties for 
nonmagnetic or collinearly magnetized electrodes. When 
dealing with noncollinear magnetizations, a new physical 
effect comes into play: the precession of the dot spin in 
an exchange field generated by virtual tunneling between 
the dot and the leads. As mentioned in the introduction, 
the interplay between spin accumulation and spin pre- 
cession gives rise to a number of characteristic transport 
signatures. In the following, we discuss how the depen- 
dence of the exchange field on the state of the impurity 
spin influences the transport properties, in particular, the 
current and the finite-frequency Fano factor. 

1 . Current 

In the normal quantum-dot spin valve, the current- 
voltage characteristics provides information about the ex- 
change field acting on the dot spin. At the particle-hole 
symmetric point of the Anderson model (e = —U/2), the 
exchange field vanishes. Hence, only the spin accumula- 
tion on the dot influences the current, leading to a re- 
duced current in the noncollinear geometry compared to 
the parallel case. Away from this point, a finite exchange 
field acts on the dot spin and gives rise to a precession of 
the spin out of the blocking position. As a consequence, 
the current flowing through the system becomes larger 
than at the particle- hole symmetric point, resulting in a 
characteristic U-shaped I — V curve with a broad region 
of negative differential conductance. 

For the quantum-dot spin valve with an impurity em- 
bedded in one of the tunnel barriers studying the current 
as a function of bias voltage allows us to investigate the 
dependence of the exchange field on the state of the im- 
purity. We will first focus on the situation where the 
interference term has positive sign and V > 0. In this 
case, if the impurity is in state | + S), there is a large 
exchange field acting on the dot spin (cf. Fig. 3), giving 
rise to a strong precession and thereby to a clear lifting 
of the spin blockade. If the impurity is in state | — S), 
however, there is only a small exchange field acting on 
the dot. Consequently, the spin blockade persists and 
the current is not much enhanced by the exchange field. 

As is shown in Fig. 9(a) (solid black curve), we there- 
fore find a large current at the onset of transport through 
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FIG. 9. (Color online) Current as a function of bias voltage 
for perpendicular magnetizations for p = 0.9 and IY = Tr = 
10 Jr. Black curves are for 77 = +1 while red curves are for 
77 = —1. For the upper panel, e = — 300fcB-T while for the 
lower panel e = 50/cbT. Other parameters as in Fig. 3. While 
the solid lines represent the full result, for the dashed curves 
the exchange field was set to zero by hand. 



the quantum dot. When increasing the bias voltage, 
the current decreases slightly as the exchange field be- 
comes weaker. When the bias voltage is increased above 
the threshold for impurity excitations, the impurity is 
switched from its ground state | + S) into the state | — S) 
by the spin-polarized current (cf. upper panel of Fig. 10). 
As a consequence, the exchange field changes abruptly at 
threshold towards smaller values. This implies that the 
spin blockade on the dot cannot be lifted anymore and 
the current is suppressed compared to its values below 
threshold. 

When the interference term has a negative sign [red 
curves in Fig. 9(a)], the situation is reversed. Now the 
exchange field takes on small values when the impurity 
is in state | + S) while it takes large values when the 
impurity is in state | — S) . Hence, the current is now 
suppressed below threshold while above threshold we find 
a nontrivial bias dependence of the current due to the 
energy dependence of the exchange field. 



0.4 
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FIG. 10. (Color online) x, y and z component of the average 
dot spin and 2 component of the average impurity spin for 
perpendicular magnetizations as function of bias voltage for 
positive (upper panel) and negative (lower panel) sign of the 
interference term. Parameters as in Fig. 9(a). 



When a negative bias voltage is applied, no switching 
occurs at the impurity excitation threshold as the spin- 
polarized current has the tendency to bring the impurity 
into the state | + S) which is the ground state. We there- 
fore find that now the exchange field effects are clearly 
visible for a positive sign of the interference term, while 
they are very small for a negative interference term. 

In summary, we have shown how the current-voltage 
characteristics for noncollinear magnetizations provide 
access to the dependence of the exchange field on the 
impurity spin state as well as to the sign of the inter- 
ference term. We note, however, that the observation of 
exchange field effects in the current relies on large polar- 
izations. 



2. Frequency- dependent Fano factor 

The effects of the exchange field switching on the I — V 
characteristic occur only for large polarizations as oth- 
erwise the spin blockade on the dot is too weak. An- 
other method to gain information about the exchange 
field which also works for smaller polarizations is to study 
the frequency-dependent Fano factor. 35 

In Fig. 11, the frequency-dependent Fano factor is 
shown for different bias voltages and different parameters 
of the spin Hamiltonian. Parameters are chosen such that 
for the smaller of the two bias voltages (black curves) the 
impurity spin cannot be excited. For the larger bias volt- 



16 




V 


= 7()k B T, B 


= 25fc B T 


V 


= 80fc B T, B 


= 25fc B T 


, V 


= 70fc B T, B 


= 50fc B T 


, V 


= 80fc B r, B 


= 50fc B T 



p = 0.3 





FIG. 12. (Color online) Schematic of the conductance for the 
situation where B » T, J < T. The lines indicate peaks in 
the differential conductance. 
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FIG. 11. (Color online) Frequency-dependent Fano factor for 
moderate (p = 0.3, upper panel) and large (p = 0.9, lower 
panel) polarizations. For p = 0.3, an external magnetic field 
perpendicular B cxt = 1.5F is applied to the quantum dot 
alongs its quantization axis. Other parameters are the same 
as in Fig. 9, except that now D = and rj — +1. 



age, the impurity parameters allow an excitation of the 
impurity spin only for the red solid curves while for the 
red dashed curves the impurity still stays in the ground 
state. 

For large polarizations and voltages below the excita- 
tion threshold, the finite-frequency noise shows a peak at 
the Larmor frequency of the exchange field. As the im- 
purity is in the ground state, the system is not sensitive 
to the impurity parameters and both curves practically 
coincide. If the bias is increased but the impurity stays in 
the ground state, the absolute value of the exchange field 
is slightly reduced, cf. lower panel of Fig. 3, resulting in 
a small shift of the resonance signal toward smaller fre- 
quencies. If the bias is increased and a switching of the 
impurity spin occurs, we instead find that the exchange 
field is significantly reduced, cf. Fig. 3, and therefore the 
resonance peak is also shifted to much smaller frequen- 
cies. Hence, by detecting the resonance frequency as a 
function of bias voltage one can gain information about 
the switching of the exchange field as a consequnce of the 
switching of the impurity spin. 

For small polarizations, the effect does not work the 
same way as just described. Now, the exchange field just 
changes sign upon switching the impurity spin, cf. Fig. 3. 
Therefore, no clear shifting of the resonance position oc- 
curs. To overcome this problem, an external magnetic 
field of strength comparable to the exchange field can be 
applied perpendicularly to the plane defined by the elec- 



trode magnetizations. Now, the position is defined by 
the Larmor frequency of the total magnetic field while 
the form of the resonance signal depends on the relative 
angle between external field and exchange field as can be 
seen in the upper panel of Fig. 11. If the impurity spin 
stays in the ground state, a peak shows up in F(lu) while 
a shoulder occurs at the Larmor frequency if the impu- 
rity spin can be switched by the current. Hence, we have 
shown that the switching of the exchange field can be 
monitored also for moderate polarizations by measuring 
the finite-frequency current noise for a series of applied 
bias voltages. 



VI. RESULTS - SMALL SPIN ON THE DOT 

In this section, we discuss the transport properties of a 
quantum-dot spin valve containing an additional spin 1/2 
on the quantum dot. We will focus on the two transport 
regimes (ii) and (v) introduced above, cf. Table II. In 
these regimes, the exchange coupling is small, J < T, 
making them particularly suited to describe the influence 
of nuclear spins on transport through a quantum-dot spin 
valve. In regime (ii) , the externally applied magnetic field 
B is much larger than the tunnel coupling between dot 
and leads. In this regime, we show that the coherent 
superpositions of the singlet and one of the triplet states 
do not show up in the current. However, they do lead to 
a signature in the finite-frequency Fano factor. In regime 
(v), the external magnetic field is weak, B < T. In this 
case, we discuss how to extract the exchange coupling and 
external field from measurements of the finite-frequency 
Fano factor. 



A. Large magnetic field 

We start our discussion with the case B ^> T, J < T, 
where only superpositions of the singlet S and the triplet 
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state T° are relevant. In Fig. 12 we show schematically 
the differential conductance as a function of bias volt- 
age V and level position e. The thick black lines indi- 
cate where a dot level is in resonance with either the 
left or right Fermi energy. Hence, they separate the 
Coulomb-blockaded regions from the regions where se- 
quential through the dot is possible. For concreteness, 
we will now discuss the sequence of transport processes 
that come into play upon increasing the bias voltage for 
a fixed level position e = B/2. For small bias voltages, 
transport is possible only through the states |0 1) and 
|T~) as transitions from |0 I) to any of the other singly 
occupied states are energetically forbidden while transi- 
tions from \T~) to |0 f) are impossible as they violate 
spin conservation. 

When the bias voltage is increased across the blue 
dashed line, the bias window becomes large enough to 
also allow transitions from |0 \.) to \T°) and \S). As for 
both of these states there is a finite probability to find the 
impurity spin in state f, both, \T°) and \S), can serve as 
a starting point for a transition into the state |0 t) which 
then allows transitions into state \T + ). Hence, we find 
that as soon as the bias is large enough to excite the 
|0 i)-\T°) and |0 transition, all empty and singly 

occupied dot states contribute to transport. 

Upon increasing the bias voltage further across the red 
dashed line, another set of transport processes becomes 
energetically possible. While transitions from singly oc- 
cupied dot state with the lowest energy, \T~), to the 
doubly occupied state \d J.) are still not possible, it is 
nevertheless possible to occupy the dot with two elec- 
trons by taking as a starting point either \T°) or \S). 

Finally, when crossing another black line, the bias volt- 
age is large enough to allow transitions between any two 
dot states that conserve spin. 

Approximating the Fermi function as step functions, 
one can derive the following analytical expressions for the 
current through the quantum dot in the different trans- 
port regimes: 



The Fano factor at zero frequency is given by 
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It is interesting to note that these expressions are pre- 
cisely the same as found by Thielmann et al. 86 for trans- 
port through a single-level quantum dot subject to a large 
magnetic field coupled to normal leads. Hence, in the 
chosen parameter regime, the current through the system 
and the zero-frequency Fano factor are neither sensitive 
to the presence of the impurity spin nor to the presence 
of ferromagnetic leads. The absence of any magnetore- 
sistance is due to the presence of the large external field 
perpendicular to the plane spanned by the magnetiza- 
tions which renders the system insensitive to the relative 
orientation of the magnetizations in this plane. 

The fact that the system is insensitive to the presence 
of the impurity spin and to the coherent superpositions 
between \S) and \T°) deserves some further investiga- 
tion. Typically in systems where coherent superpositions 
of different states have to be taken into account as, e.g., a 
normal quantum-dot spin valve 32 ' 33 or double quantum 
dots, 8 the current shows a nontrivial bias dependence 
with broad regions of negative differential conductance 
due to the energy-dependent level renormalization that 
arises from virtual tunneling between the dot and the 
leads. Such effects are clearly absent here. To under- 
stand this behavior, we analyze the term describing the 
isospin accumulation, Eq. (4.9), in more detail. Due to 
the Fermi functions in this expression, a finite isospin ac- 
cumulation can arise only in the regions marked III in 
Fig. 12. When plugging in the expressions for the occu- 
pation probabilities, we find however that 
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L 4(r L + r R ) 2 (r 2 + 2r L r R + 2r 2 ) 
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i.e., the contributions to the isospin accumulation from 
the left and right lead cancel precisely. In consequence, 
the isospin vanishes on average in the stationary state 
such that the aforementioned level renormalization ef- 
fects cannot influence the current and zero-frequency 
noise. 

So far, it seems that the coherent superpositions of 
IS*) and |T°) do not influence the transport through the 
system at all, such that a description in terms of ordi- 
nary rate equations which neglects the coherences suf- 
fices. This is not true, however, as can be seen from the 
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FIG. 13. (Color online) Finite-frequency Fano factor for dif- 
ferent values of the exchange coupling J. The Fano factor 
shows a peak at the Larmor frequency associated with the 
exchange field. Parameters are Fl = 2Fr, fc B T = IOFl, 
B = 50fc B T, U = 5B, e = B/2 and V = 



finite-frequency Fano factor, which is shown in Fig. 13. 
Here, apart from the usual peak which arises at zero fre- 
quency, the Fano factor additionally shows a peak asso- 
ciated with the Larmor frequency of the exchange field in 
Eq. (4.11). The coherent superpositions play a role here 
as the finite-frequency noise is sensitive to the dynamics 
of the dot spins while the current only captures stationary 
properties of the quantum dot. While the contributions 
form the left and right lead to the isospin accumulation 
cancel on average, indivdual processes can nevertheless 
give a finite isospin accumulation. 

In contrast to the normal quantum-dot spin valve, 
where the Larmor frequency associated with the ex- 
change field depends on the level position and applied 
bias voltage, 35 we find that in the system under inves- 
tigation here, the Larmor frequency is simply given by 
the exchange coupling strength J. This is due to the fact 
that the energy-dependent contributions to the exchange 
field arise only in the x component which is parallel to the 
isospin accumulating on the dot and hence cannot influ- 
ence its precessional motion. The frequency-dependent 
Fano factor may, therefore, be used to experimentally 
determine the exchange coupling. 



B. Small magnetic field 

In the last section, we saw that coherent superposi- 
tions can give rise to a signal in the frequency-dependent 
Fano factor F(u) which could be interpreted as the pre- 
cession of the isospin with the Larmor frequency of the 
exchange field. In the case where B,J<T, the situation 
is more complicated as now superpositions between any 
two states with the same number of electrons on the dot 
have to be taken into account. 

In Fig. 14, we show the finite-frequency Fano factor 
F(uj) as a function of the frequency u and the externally 
applied magnetic field B for fixed exchange coupling J 




FIG. 14. (Color online) Frequency-dependent Fano factor as 
a function of frequency and external magnetic field (upper 
panel) or exchange coupling (lower panel) . The level splittings 
calculated for a system of two exchange coupled spins in a 
magnetic field which is the sum of the exchange fields and the 
external fields are indicated by dotted lines. Parameters are 
J = 2T (upper panel), B = 2F (lower panel), T L = 10r R = 
0.05fc B T, e = 50fc B T, V = 105fc B T, U = 150fc B T, p = 0.9, 
= tt/2. 



and vice versa. In both cases, the Fano factor shows a 
number of features at frequencies which all show a non- 
trivial dependence on the external field B and the ex- 
change coupling J. To gain a better understanding of 
these features, we consider a simpel spin model for the 
quantum dot. This is reasonable as the dot is singly oc- 
cupied most of the time for the chosen parameters due 
to the asymmetric coupling to the leads. We model the 
dot as consisting of two spin 1/2 particles that are ex- 
change coupled and subject to external magnetic fields. 
While the impurity spin only couples to the externally 
applied field, the electron spin additionally experiences 
the exchange field generated by quantum charge fluctu- 
ations on the dot. Hence, the Hamiltonian of our model 
is given by 



H — JSi • S2 + Bi ■ Si + B2 ■ S2 



(6.10) 



where Bi = B CXj LnL + -B C x.R.nR + Be z is the sum of 
the external magnetic field and the exchange field while 
B2 = Be z is simply given by the external field. To show 
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that this is indeed the correct way to describe the quan- 
tum dot system, in Appendix D we give the equation of 
motion for the density matrix elements of the two spins. 
Comparing them to the master equation for the quantum 
dot system, Eqs. (4.14)-(4.19), shows the correspondence 
between the two systems. The only difference is the ab- 
sence of dissipative terms in the equations for the spin 
model. 

The differences between the energies of the various 
eigenstates of Eq. (6.10) as a function of exchange cou- 
pling and external field are indicated as black dotted lines 
in Fig. 14. We find a nice agreement between these energy 
differences and the features observed in the Fano factor. 
On the one hand, this indicates that the finite-frequency 
Fano factor is indeed sensitive to level splittings compa- 
rable to the tunnel couplings which cannot be resolved 
in the differential conductance as the conductance peaks 
are broadened by temperature which satisfies T <C fceT. 
This can provide experimental access, e.g., to the cou- 
pling between the dot and impurity spin. 

On the other hand, our results nicely demonstrate that 
the exchange field which arises due to virtual tunneling 
processes and not due to stray fields from the ferromag- 
netic leads, acts only on the electron spin but not on the 
impurity spin. If one could measure the Fano factor at 
finite frequency as computed above, one could therefore 
clearly distinguish exchange field effects from stray field 
effects as the latter ones would act in the same way on 
both spins and thereby give rise to a trivial dependence 
of the energy differences on the exchange coupling and 
the external field. 



VII. CONCLUSIONS 

In this paper, we investigated transport through 
quantum-dot spin valves containing magnetic impurities. 
First, we considered the case where a large, anisotropic 
impurity spin is localized in one of the tunneling bar- 
rier. Our main focus was on the spin excitation and spin 
switching on the impurity spin. We pointed out how 
the Coulomb interaction on the dot allows a more de- 
tailed spectroscopy of the impurity spin. For magnetic 
electrodes, the spin-polarized current through the sys- 
tem can switch the spin of the impurity. As the state 
of the impurity spin influences the current, an interplay 
between the current and the impurity arises which gives 
rise to a series of positive and negative differential con- 
ductance as a function of bias voltage. We, furthermore, 
found that the dynamics of the impurity spin gives rise 
to a random telegraph signal associated with a huge Fano 
factor. Additionally, we found that for tunneling through 
the barrier containing the impurity interference between 
direct and exchange tunneling takes place even for non- 
magnetic systems. This is in contrast to tunnel barriers 
where such interference effects arise only for magnetic 
electrodes. Finally, we found that the exchange field act- 
ing on the dot spin becomes dependent on the impurity 



spin state. This dependence allows the observation of 
the impurity switching in the current as well as in the 
finite-frequency noise for noncollinear magnetizations. 

Second, we considered a spin 1/2 impurity localized on 
the dot in the regime of small exchange coupling to the 
electron spin. Here, our main focus was on the coherent 
spin dynamics on the quantum dot. For a large exter- 
nal magnetic field, only coherent superpositions between 
the singlet \S) and triplet |T°) are relevant. However, 
they influence neither the current nor the zero-frequency 
noise as they vanish in the stationary state. Neverthe- 
less, as the superpositions can be excited, they give rise 
to a signal in the finite-frequency noise at u> = J. For 
small external magnetic fields, all coherences have to be 
taken into account. They give rise to a large number of 
features in the finite-frequency noise. These provide in- 
formations about the level splittings that are influenced 
by renormalization effects due to virtual tunneling to the 
leads. 

To summarize, we found that in both cases, for the 
large as well as for the small impurity spin, there is an in- 
tricate interplay between the electron and impurity spin 
that manifests itself in the transport properties of the sys- 
tem. The details of of this interplay vary, however. The 
large impurity spin does not show any coherent dynamics 
in the sequential tunneling regime as its level splittings 
generated by the anisotropy are large compared to the 
tunnel couplings and the temperature. In consequence, 
these splittings can be resolved in the current which also 
allows to study the current-induced switching of the im- 
purity. In contrast, for a small impurity spin, both, the 
electron and the impurity spin behave as one object that 
displays a coherent dynamics as level splittings are of the 
order of the tunnel coupling. Clear signatures of this dy- 
namics can be found in the finite-frequency noise that 
allows to resolve level splittings of the order T <C k^T. 
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Appendix A: Diagrammatic rules 

In this appendix, we summarize the diagrammatic 
rules to evaluate the self-energies entering the general- 
ized master Eq. (3.1) and the expressions for the current 
in Eq. (3.4) and current noise in Eq. (3.5). To compute 

the self-energy W XlX } , 
X2X2 

1. Draw all topological different diagrams with tun- 
neling lines connecting vertices on either the same 
or opposite propagators. Assign to the four corners 
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and all propagators states x an d corresponding en- 
ergies E x as well as an energy uj for every tunneling 
line. 

2. For each part of the diagram between adjacent ver- 
tices, assign a resolvent 1/(AE + i0 + ) where AE 
is the difference between the energies of left- and 
right-going tunneling lines and propagators. 

3. For each tunneling line, the diagram acquires a fac- 
tor of J^f^^ui), where the sign is determined by 
whether the line runs forward (-) or backward (+) 
with respect to the Keldysh contour. 

4. For each pair of vertices connected by 
a tunneling line the diagram is multi- 
plied by l+P pr ( x [\C rt \ Xi )( x ' f \ci t \ X{ ) + 

¥Pr(x!i\Cri\xi)(x!t\Cti\Xi), where Xi and X [ 
(Xf and x'[) are the states that enter and leave 
the vertex the tunneling line begins (ends) at, 
respectively. The operators C r ^ and C r ± are the 
coefficients (including the dot operators!) of a^ k+ 

and al k _ in the tunnel Hamiltonians (2.5), (2.6) 
and (2.9), respectively. 

5. Each vertex connecting the states d and 4- gives rise 
to a factor of —1. 

6. The diagram obtains a factor (— l) a+b , where 



a is the number of vertices on the lower propagator 
and b is the number of crossings of tunneling lines. 

7. Sum over all leads r and integrate over all energies 

LO. 



ixix'i 
xix'i 



and W X1X , 1 are obtained by 

X2X2 J 



The diagrams for W 
replacing one, respectively two tunneling vertices by cur- 
rent vertices. These give rise to a factor of +1/2 (-1/2) 
if the vertex is on the upper (lower) branch and corre- 
sponds to an electron tunneling into the right (left) lead 
or out off the left (right) lead. 



Appendix B: Master equation 

1. Model A: Large spin in the barrier 

In this appendix, we give explicit expressions for the 
various quantities that occurred in the master equations 
for the occupation probabilities and the spin, Eqs. (4.2) 
and (4.3) of model A, i.e., a quantum-dot spin valve with 
an impurity spin embedded in the right tunnel barrier. 
Introducing 7 = Fr + m 2 J R + 2mpr)^/T R J K , 7 = pT R + 
m 2 pjR+2TO?7-v/r R J R , a± = B+(2m±l)D, and A±(m) — 
S(S + 1) — m(m ± 1), we find for the coupling to the left 
lead 
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which is identical to the expressions found for the ordinary quantum-dot spin valve in Ref. 33. For the coupling to 
the right lead, we have 



( <oo 7/r(£) 

27/r» Wi% 2 7 / R " (e + U) 
V lf + {e + U) W' 



(B5) 



(0) 
R,dd 



where the diagonal matrix elements are determined by the sum rule e T W = 0, 



Vi° ] = 2 



( \ 

7 + /r (e + E0] + ^+(™ - [^/r (e - «-) + V/r ( £ + 17 + «-)] 

-A_(m + 1) J R [ V/r (e + «+) + ^/r (e + ^ - «+]) 
V -7/+(e + C7) / 



i±E/- (e _ a+ ) 

W R +1) =A_(m + l)J R | (l-p)/+( e + a+) (1 + p)/ R (e + C/ - a+) 

V/R( £ + (7 + a +) 



V^ +1) =A_(m + l)J R 



-(l+p)/-( £ -a+) \ 


-(l-p)/+(e + tf + a + )/ 



V/R( £ + a -) \ 

W R T 1) =A + (m-l)J R I (l+p)/+( e -a_) " (1 - p)/ R ( £ + U + a_) 

l+£/+( £ + c/_ a _) o y 



V R 1} = A+(m-l)J R | () 



(l-p)/K(e + a_) \ 


{l+p)f+{e + U-a-)J 



(o) 

acc,R 



= \ WOO^o,™ - 7/r (e + ^)^,m + 

Jr /1+P 



2 

- 1 > < - 1)^ ( £ - «-) + ( £ + + «-) 

-A_(m + 1)^ (^/r (e + «+) + ^/r"^ + ^ - «+) 



L,m | 



-Pl.m f n R 



(+i) 



acc.R 



A_(m + 1)J R 



(1 - p)f+(e + a+)P , m +i + (1 +p)/r (e + - a+)P d , m+ i 



n R , 



"(ft" 



(-1) 



acc.R 



-A+(ra- 1)J R 



(1 - a-)P , m -i + (1 - p)/r (e + £/ + a_)P d , m _i 



n R , 





r - 




' rel.R 



7(/ R (e) + Jr~ (e + E0) + ^+(™ - 1) Jr. (^/r 0* - «-) + ^/r (* + ^ + «-) 
+.4_(m + 1) J R (^/r (e + a+) + ^/r + ^ - «+) 
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2. Model B: Small spin on the dot 



The vector V occurring in the master Eq. (4.7) of model B, i.e., the quantum-dot spin valve containing a spin 1/2 
on the dot is given by 

/ -f-(e-B/2) \ 

f-(e + B/2) 


f-( £ -B/2)-f-(s + B/2) + f+(e + U + B/2)-f+( S + U-B/2) . (B15) 



-f+(e + U + B/2) 
V f+(e+U-B/2) J 



v = E r - 



r 



Appendix C: Density matrix elements 

In this section, we give the relation between the phys- 
ical quantities introduced in Sec. IV B 2 and the density 
matrix elements in Eq. (4.6). For the occupation proba- 
bilities of the dot we have 

P = P at + Pou (CI) 
Ft = P T+ + p T0 + p T _ + p Sj (C2) 

Pd = Pd\ + Pdi- (C3) 

For the expectation values of the impurity spin when the 
dot is empty, we get 



The expectation value of the electron spin on the dot can 
be expressed as 



Six 

Sly 

Slz 



Re P§_ - Re P$ + + R e P%_ + Re P% + 



V2 



(C13) 



Im P|_ - Im P§ + + Im P^l - Im P*"! 

! (C14) 



V2 

P T + -Pt- +2ReP^ 



(C15) 



For the scalar product of the electron and impurity spin, 
we find 



Si • S 2 — 



P T + + P T + P 7 



3Ps 



(C16) 



Sox 


— Re Pq^ , 


(C4) 


Soy 


= ImP°|, 


(C5) 


Sqz 


_ Pot - Poi 
2 


(C6) 



When the dot is singly occupied, we have 

- Re P£_ + Re P§ + + Re P^! + Re Pt 



S 2 



Im P§_ + Im P| + + Im P^l - Im P^! 

S 2 y 
S 2Z 



T+ 



V2 

P T+ -P T - -2ReP| 



and for the doubly occupied dot, we find 



Sdx — RePfi, 



(C7) 
(C8) 
(C9) 



Sdy 



S, 



dz 



Tm P d ^ 
1111 r d\ ' 

Pdt — Pdj 

2 



(CIO) 
(Cll) 

(C12) 



(Si x S 



2)y 



V2 

ReP£ + + ReP|_ 



V2 



(Si x S 2 ) z = ImP# . 
Finally, the quadrupole moments are given by 
2P T o - P T - - P T + + 6ReP£+ 



Qxx ~ 
Qyy 

Q zz 
Qxy 
Qxz 
Qyz 



12 

2P T o - P T - - P T + - 6RePf+ 
12 

P T + - 2P T o + P T - 



6 



ImP^ + 



ReP^! + RePf° 



2^2 
+ Ii 
2\/2 



ImP£_ +ImP^ + 



(C17) 

(C18) 
(C19) 



(C20) 

(C21) 
(C22) 

(C23) 

(C24) 

(C25) 
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Appendix D: Equation of motion for two spin 1/2 
particles 

The equations of motion of two exchange-coupled spins 
Si and S2 subject to magnetic fields Bi and B2, respec- 
tively, as described by the Hamiltonian (6.10) are given 
by ' 

^ = -SixBi-J(SixS a ), (Dl) 



— ^ = -S 2 xB 2 + J(S 1 xS 2 ), (D2) 



^(S 1 -S 2 ) = (S 1 xS 2 )-(B 1 -B 2 ), (D3) 



^(Si x S 2 ) = -1(S X x S 2 ) x (Bi + B 2 ) + (q - |(Si • S 2 )) • (B x - B 2 ) + ^(S x - S 2 ), (D4) 



d 
~dt 



1 /1 



2 v 2 



(Si x Sa) 4 (B x - B a )j + J(S X x S 2 )i(B x - B 2 ) 4 - i(Si x S 2 ) • (B : - B 2 )<J i; 



:^iimQij(Bi + B 2 ) m — 7r£jimQii(Bi + B 2 ) 



(D5) 



It is interesting to note that some terms contain the sum 
of the two magnetic field while other couple to the differ- 
ence. In consequence, for the quantum-dot system under 
investigation, some terms are only sensitive to the ex- 



change field as the external field acts on both spins in 
the same way, while other terms are sensitive to both 
field contributions. 
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